Resonant amplification of the Andreev process in ballistic Josephson junctions 
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We study the Josephson effect in balhstic SfNfS and SfFfS double-barrier junctions, consisting 
of superconductors (S), a clean normal (N) or ferromagnetic (F) metal, and insulating interfaces 
(I). For short SINIS double-tunnel junctions with one channel open for quasiparticle propagation, 
the critical Josephson current as a function of the junction width shows sharp peaks because of a 
resonant amplification of the Andreev process: when the quasibound states of the normal interlayer 
enter the superconducting gap the Andreev bound states are lowered down to the Fermi level. For 
corresponding SIFIS junctions the quasibound states are spin-split; they amplify the supercurrent 
less efficiently, and trigger transitions between and tt states of the junction. In contrast to SI- 
NIS junctions, a narrow dip related to the - tt transition opens up exactly at the peak due to a 
compensation of partial currents flowing in opposite directions. With an increased barrier trans- 
parency these features are gradually lost, due to the broadening and overlap of quasibound states. 
Temperature-induced transitions both from to tt and from tt to states are studied by computing 
the phase diagram (with temperature and junction width as the variables) for different interfacial 
transparencies varying from metallic to the tunnel limit. 
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PACS numbers: 74.50.-|-r, 74.45.-l-c 



I. INTRODUCTION 

Recently, many interesting phenomena were investi- 
gated in superconductor (S) - ferromagnet (F) Josephson 
junctionsiSi One of the most interesting effects is the so- 
called TT state of SFS junctions, in which the ground state 
is characterized by an intrinsic phase difference (f) = tt 
between two superconductorsi^i^i^iSiL2i2ii£ The tt junc- 
tions could also have important applications, such as for 
quantum qubitsiii The Josephson effect in double-barrier 
SIFIS junctions, with a clean metallic ferromagnet be- 
tween two insulating interfaces (I), has been studied by 
solving the scattering problem based on the Bogoliubov 
- de Gennes equation; it was found that spin-split quasi- 
particle transmission resonances contribute significantly 
to the Andreev process^^^i^ A comprehensive numerical 
study of clean layered S /F heterostructures has also been 
given recently^ 

For SINIS structures (with a clean non ferromag- 
netic normal metal N), sharp peaks occur in the criti- 
cal Josephson current Ic as a function of the junction 
width d for each transverse channel. This is the result 
of a resonant amplification of the Andreev process when 
the quasibound states in the INI interlayer approach the 
Fermi level. This effect is the superconducting analog 
of the Bohm resonant tunneling,M previously studied 
in RefsiiSiiLiSiiS For corresponding SIFIS junctions, the 
resonant nature of - tt transitions comes into play in 
the tunnel limit. Triggering of the - tt transitions by 
resonant amplification, as well as the coexistence of sta- 
ble and metastable and tt states in the vicinity of the 
transition, can be fully understood. In contrast to 
SINIS junctions where the critical current reaches a peak 
value when the Andreev bound states cross the Fermi 



level, here a narrow dip related to the - tt transition 
opens up exactly at the peak due to the compensation 
of partial currents flowing in opposite directions. We 
emphasize that this is not the case for a higher interface 
transparency due to the broadening and overlap of quasi- 
bound states, nor for planar junctionsJS Experimentally, 
SIFIS nanostructures could be realized, for example, in 
gated heterostructures,^^ ferromagnetic nanoparticles^^ 
or in the break junctionsSi in external magnetic field pro- 
ducing the Zeeman splitting of Andreev levelsiS^ 



It has become a paradigm that in the tt state of SFiS 
junctions, with a thin barrier of insulating ferromagnet 
(Fi), there is net spin polarization of current-carrying 
Andreev levels However this "rule" is not quite 
general. In this paper we show that in low-transparent 
SIFIS nanostructures with one open transverse channel 
the situation can be the opposite. We can generalize the 
previous rule: the - tt transition is accompanied by the 
switching of the net spin polarization of current-carrying 
Andreev levels. The critical current oscillates as a func- 
tion of the junction width, and the resonant amplifica- 
tion occurs when quasibound states of the ferromagnet 
between two insulating layers enter the gap. The - 
TT transitions take place when Andreev levels intersect 
at the Fermi energy, being accompanied by a switching 
of their spin polarization. Experiments show that the 
TT state of a dirty SFS junction can correspond to lower 
temperatures than the state. ^ The theory of ballistic 
one-barrier SFS junctions predicts the opposite We 
demonstrate that in the phase diagram of double-barrier 
ballistic SIFIS junctions the tt state can correspond to 
lower temperatures than the state, depending on pa- 
rameters. 
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II. JOSEPHSON CURRENT 

We consider a ballistic double-barrier Josephson junc- 
tion: a normal-metal (N) or ferromagnetic (F) interlayer 
of thickness d is connected to superconductors (S) by in- 
sulating non-magnetic interfaces (Ii, I2), see Fig. ^ We 
assume that both metals are clean and that the left and 
the right superconductors are equal with the constant 
pair potential Aexp(±z0/2), where (j) is the macroscopic 
phase difference across the junction. Within the Stoner 
model, the ferromagnet is characterized by the exchange 
potential h, and the insulating interfaces are modeled by 
(5-function potential barriers 

The Josephson current can be expressed in terms of 
the free energy of the junction as 



where 7 is defined by 

cos 7 = r+r_ cos /? + cos (/). (5) 

Here, a = arccos(i?/A), ^ — x+ — X-, and (3 = (x+ — 
X+) — (X- ~ X-)- From Eq. Q the Andreev levels are 

^^7,^(0) = Acos ( IsignsmI 1 , 

(6) 

where $0 = = 0,a) and 70 = j{E = 0,a). 

In the limiting case of a single insulating barrier this 
formula reduces to the corresponding equations found 
in Refs. p7,28^29j^ Using Eq. |5j in the form d^j = 
t+f_ sin(/)/sin7, we finally obtain 
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where the phase-dependent part of the junction free en- 
ergy, F{(j)), is given by 



2cosh 



2kBT 



(2) 



Here, Ena{<P) are the Andreev levels measured with re- 
spect to the chemical potential, n denotes a transverse 
channel, and a = ±1 is the spin quantum number. 

Assuming the SIFIS junction shorter than the su- 
perconducting coherence length, we can calculate the 
Josephson current using the relation 
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where fna{T) is the Fermi distribution, and En„{(j)) are 
zeroes of the common denominator of the scattering am- 
plitudes G , calculated from the Bogoliubov - de Gennes 
equation and given by Eq. (17) in Refii^ 

In the following, we assume that only one channel 
(n 1) is open in the short superconducting junction. 
Without a ferromagnetic region there are two spin degen- 
erate Andreev levels with the same absolute value of en- 
ergy, but with the opposite sign that can be parametrized 
by rj = ±1. The ferromagnet spin-splits Andreev levels, 
and there are four of them, Eria{4>)- We will express 
Ejj„{(j)) explicitly in terms of the scattering matrix. ^^-S^ 
For the IFI structure formed by the ferromagnet and two 
insulating barriers we define t± and r± to be the absolute 
values of the electron and hole transmission and reflec- 
tion amplitudes, t{ztE,a) and f{zLE,a), where E is the 
quasiparticle energy with respect to the chemical poten- 
tial. The corresponding phases are denoted as x± and 
X± ■ The Andreev levels can be calculated as roots of the 
equation 



cos(<I> — 2a) — cos 7 = 0, 



y eAM-Bin0^. / 70 + ^^0 
^^^^^2h sin7o I 2 



2kBT 
(7) 

which is a generalization of Eq. (3) in Ref. fSi] . The trans- 
mission amplitudes of the IFI structure can be expressed 
through the transmission and reflection amplitudes of the 
insulating barriers, i/j^ ^ , f/^ ^ , and the momentum of an 



electron in F, qa{E) 
k = 



2m{EF 



E + (j/i), and in S, 



2m{EF'-'^^ + n), where Q = ^fW~~^. Namely, 



t{E,G) 

f{E,a) 



til tl2 ^ 



i{q„—k)d 



1 — fj^f/^e^*^'''' ' 



-ikd 



1 — r'j_^ri^e^'^^''^ 



(8) 



Note that the expression for the transmission amplitude 
has the same form as the one for the Fabry-Perot interfer- 
ometer. Here, the scattering amplitudes of the insulating 
barriers are 



til,: 



k' 
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q„ + 2iZ^,2kP ' 

- q„ - 2iZi^2k^P 

- q„ + 2iZi.2fc(/' ' 



(9) 



where k 



(S) 



2mEp^ the parameters Z12 deter- 
mine transmission probabilities of insulating barriers, 
and r' = —r*t/t* (* denotes the complex conjugation). 
In the tunnel limit, l^/j 2 1 ~* 0, the positions of the bound 
states are simply given by the conditions dqa-{+E) = niir 
and dqa{—E) = n27r, where ni and n2 are integers. 

The results are illustrated in Figs.EHTTl In all illustra- 
tions superconductors are characterized by A/Ep^^'' = 
10"'^, the Fermi energies of the two metals are equal, 
Ep^^\ in F or N interlayers the normalized 

0.1 or 0, and the 



E, 



(4) exchange energy is A" = h/Ep 
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insulating barriers are assumed to be equal and char- 
acterized by the common parameter = — 
related to the one-barrier transmission probability via 

In SINIS double-tunnel junctions sharp peaks occur in 
the critical Josephson current as the junction width is 
varied (Fig. [3 top panel) because of resonant amplifica- 
tion of the Andreev process when the quasibound states 
of the INI interlayer enter the superconducting gap and 
lower the phase-sensitive Andreev bound states down to 
the Fermi level (Fig. O . The order of magnitude of the 
peaks is the same as the critical Josephson current val- 
ues in transparent junctions. Between two peaks, the 
junction behaves as a usual SIS junction. In this case, 
quasibound states are not in the superconducting gap. 

In SIFIS junctions the critical current significantly de- 
creases if the junction transparency is reduced (Fig. 12 
bottom panel) . Instead of the critical current reaching a 
peak value, in the case of low transparency a narrow dip 
opens up exactly at the peak due to the onset of the - tt 
transition. Note that the modulation of spike heights in 
Ic (d) is due to the modulation of the Zeeman splitting. 
When the Zeeman-split quasibound states are closer, the 
amplification is larger, the peaks tend to the SNS value, 
and there is no triggering of the - tt transition. There- 
fore, only the large enough spin splitting of sharp enough 
quasibound states can trigger the - tt transition. At 
the transitions, spin-unpolarized Andreev states (with 
opposite spin orientation on the same side of the Fermi 
surface) become spin polarized, and vice versa. This is 
shown for two neighboring peaks in Fig. 21 We point 
out that, in contrast to the one-barrier problem fSL^^iSS 
or transparent interfaces, neither spin-polarized Andreev 
states have to correspond to the tt state of the junction, 
nor spin-unpolarized to the state. The following gener- 
alization is valid: every - tt transition is accompanied 
by a switching of the spin polarization of Andreev bound 
states, Fig.|3 For low transparency, switching of the spin 
polarization is not always accompanied by - tt transi- 
tions. Partial contributions of Andreev states to I{(f>) 
are shown in Fig. In the transition region, there is a 
coexistence of stable and metastable and tt states due 
to a considerable contribution of the second harmonic in 
I {(/)), Fig. EpSi In the tunnel limit, Z ^ oo, the Andreev 
bound states that are practically (j) independent coincide 
with bound states of the corresponding isolated normal- 
metal interlayer, but with nonmatching spin polarization. 
Because of this, the critical current is extremely weak. 

With increased barrier transparency, the mechanism of 
- TT transitions, described above, is modified by a broad- 
ening and overlap of quasibound states in the normal in- 
terlayer. In the case of complete transparency, the criti- 
cal current modulation is related to only two transitions 
in the considered interval of dkp, Fig. |21 bottom panel. 
This modulation is due to a finite width of quasibound 
states even for transparent interfaces because of the mis- 
match of Fermi energies for different spins resulting in a 
finite normal reflection. -° Note that also for finite, but 



relatively large transparency, Z ~ 1, there are the same 
two - TT transitions as in fully transparent junctions, 
although the critical current exhibits rapid oscillations 
due to the resonant amplification. The same result holds 
for the corresponding planar junctions after summation 
over transverse channels, even for a small transparency 
of interfacespi^ 

Temperature-induced transitions both from to tt and 
from TT to are studied by computing the phase diagram 
with temperature and junction width as variables, for 
different interfacial transparencies from metallic to tun- 
nel limit. Fig. In the case of low transparency there 
is a - TT transition for each d where the quasibound 
states cross the Fermi surface. The characteristic non- 
monotonic dependence of Ic{T) with a well-pronounced 
dip at the transition is shown in Figs. |H1 and with the 
corresponding parts of the phase diagram in the insets. 
One can see that the tt state can correspond to lower 
temperatures than the state. In Figs. 1101 and 1111 An- 
dreev levels and their contributions to /(</>) are shown for 
three temperatures. Note that the discussed mechanism 
of - TT transition is more general and robust against the 
insertion of scattering centers into the F interlayer. Sim- 
ilar results are previously obtained for a diffusive point 
contact connecting two SF bilayers (SFcFS)3 



III. CONCLUSION 

We have shown the resonant nature of - tt transi- 
tions in short ballistic SIFIS junctions with one trans- 
verse channel open for propagation. Unlike the corre- 
sponding SINIS structures, in low-transparency SIFIS 
junctions the peaks in Ic{d) curves do not reach the cor- 
responding value for transparent junctions; instead, the 
- TT transitions are triggered and narrow dips open up 
exactly at the peaks. This is not the case for higher in- 
terface transparency, nor for planar junctions. 

We have shown that the rule: "The tt state is char- 
acterized by net spin polarization of the current-carrying 
Andreev levels"^^'^® is not quite general. Remaining true 
for high-transparency junctions, in SIFIS tunnel junc- 
tions the situation can be the opposite. We generalized 
the previous rule to read as follows: - tt transition is 
accompanied by the switching of net spin polarization of 
current-carrying Andreev levels. In agreement with ex- 
perimental observations^ we also demonstrated that in 
the phase diagram of double-barrier SIFIS junctions the 
TT state can occur at lower temperatures than the state 
(or the other way round), depending on parameters. This 
also holds for planar double-barrier junctionsj^^ as well 
as for diffusive SFcFS point contactsi^ 
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Fig. 1 A schematic of a double-barrier superconductor- 
ferromagnet-superconductor junction. Ii and I2 denote 
the insulating layers. 

Fig. 2 Critical Josephson current Ic as a function of 
the normal-metal-interlayer thickness d at zero temper- 
ature for transparent (Z = 0, dotted curves) and low- 
transparency {Z = 10, solid curves) interfaces, for the 
SINIS junction {X = 0, top panel) and the SIFIS junc- 
tion {X = 0.1, bottom panel). Note that Ic{d) for low- 
transparency SIFIS junction is magnified 20 times. 

Fig. 3 A peak in Ic{d) and the spin degenerate An- 
dreev levels, E{(f>c), for a SINIS junction with Z = 10 
and X — at zero temperature. 

Fig. 4 Two neighboring peaks in Ic{d) for the SIFIS 
junction with X — 0.1 and Z = 10, and spin split An- 
dreev levels, E{(f)c)'- two spin up levels (thick and thin 
solid curves) and two spin down levels (thick and thin 
dashed curves). Insets: dips at <-^- tt transitions. 

Fig. 5 Variation of the Andreev levels as a function of 
the phase difference 4> for five values of the normal-layer 
thickness indicated in the inset of Fig. 01 with all other 
parameters being the same. The arrows indicate the spin 
polarization of the levels. 

Fig. 6 The current-phase relation (solid curves) for 
five values of the normal-layer thickness indicated in the 
inset of Fig.^ with all other parameters being the same. 
Partial currents carried by the corresponding Andreev 
levels shown in Fig.|31are also presented (dotted curves). 

Fig. 7 The phase diagram as a function of the normal- 
ized temperature T/Tc and normalized ferromagnetic- 
layer thickness dkp^ for transparent {Z — 0, top panel) 
and low-transparency interfaces {Z = 10, bottom panel), 
and for X = 0.1. 

Fig. 8 Characteristic temperature variation of Ic 
showing the temperature-induced transition from to tt 
states for a junction width dkp = 32.755, X — 0.1, and 
Z = 10. The inset shows the part of the phase diagram 
for the corresponding transition. 



Fig. 9 Characteristic temperature variation of Ic 
showing the temperature-induced transition from tt to 
states, for a junction width dkp = 59.708, X — 0.1, and 
Z = 10. The inset shows the parts of the phase diagram 
for the corresponding transition. 

Fig. 10 Variation of the Andreev levels as a function 
of the phase difference for dkp = 29.609, X = 0.1, and 
Z = 10. The arrows indicate the spin polarization of the 
levels. 

Fig. 11 The current-phase relation (solid curves) for 
the same parameters as in the Fig. ^1 and for three dif- 
ferent temperatures: (1) T/Tc = 0.01, (2) T/Tc = 0.2, 
and (3) T/Tc — 0.5. Partial currents carried by the 
corresponding Andreev levels are also presented (dotted 
curves). 
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